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Electromagnetic properties of the nucleon are explored with lattice QCD using a novel technique. 
Focusing on background electric fields, we show how the electric polarizability can be extracted 
from nucleon correlation functions. A crucial step concerns addressing contributions from the 
magnetic moment, which affects the relativistic propagation of nucleons in electric fields. By 
properly handing these contributions, we can determine both magnetic moments and electric po- 
larizabilities. Lattice results from anisotropic clover lattices are presented. Our method is not 
limited to the neutron; we show results for the proton as well. 
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1. Overview 

Electromagnetic properties of hadrons allow us to glean information about their internal struc- 
ture. Nucleon magnetic moments are well determined experimentally, and their deviation from 
Dirac's value for point-like particles, pL = Qe/2M, already implies the composite nature of the 
nucleon. The naive quark model provides the prediction: \inj\i-p = —2/3; which, while well sat- 
isfied experimentally, results from an uncontrolled approximation. Decades after the development 
of QCD as a fundamental theory of quarks and gluons, we are beginning to understand the elec- 
tromagnetism in the nucleon from first principles. A further window to the internal dynamics of 
the nucleon is given by the study of multipole polarizablities. These quantities encode the de- 
formation of the nucleon in applied electromagnetic fields, and can be accessed experimentally 
through Compton scattering. Near the chiral limit, nucleon polarizabilities arise from the deforma- 
tion of the charged pion cloud, and are highly constrained by effective interactions that emerge in 
the low-energy limit of QCD. Lattice QCD will play a crucial role in validating this low-energy 
picture, and the electromagnetic polarizabilities provide an area in which the lattice will influence 
phenomenology. 

We report on our progress in the investigation of hadronic electromagnetic properties using the 
background field method. Previously, we have shown how to determine the electric polarizabilities 
of pseudoscalar mesons using lattice QCD [1]. Recently, we have extended our analysis to the 
nucleon [2], and this work is our primary concern here. We will begin by reviewing the physics 
of spin- 1/2 particles in electric fields. In particular, we treat the spin relativistically in order to 
account for all effects at second order in the strength of the applied electric field. To determine 
the electric polarizability, we show that a background field analogue of the Born subtraction is 
necessary. We devise such a method, and apply it to the analysis of nucleon correlation functions 
calculated using lattice QCD in background electric fields. Finally we provide an outlook which 
describes refinements that must be made in order to confront experiment. 

2. Nucleon in Electric Fields 

Consider a neutron in an external electric field in Minkowski space. Because the neutron 
is composite, the external field interacts with the neutron through a tower of non-minimal cou- 
plings. For fields sufficiently weak compared to the QCD scale, \E\/Aq CD <C 1, the tower of terms 
can be ordered according to a power counting scheme. Operators with the fewest field-strength 
tensors and fewest derivatives are the most relevant in a low-energy effective field theory, while 
those with more field-strength tensors and/or more derivatives are power-law suppressed. 1 Terms 

'This theory of the neutron and photons is basically the simplest effective field theory one could imagine, however, 
there is a subtlety. In Minkowski space, there are additional non-perturbative effects that are absent in the single-neutron 
effective action. Such effects stem from the real-time production of charged pions, and arise akin to the Schwinger 
mechanism [3]. The production rate is suppressed by an exponential factor, °= exp(— jVm^/\E\), where JY is a pure 
number. This rate has been computed using chiral perturbation theory [4], but does not affect Euclidean space lattice 
simulations. Indeed, the Schwinger mechanism is absent in Euclidean space, because the essential difference between 
Euclidean electric, £ , and magnetic fields, B, is only the direction used to measure correlation functions. Consequently 
the perturbative terms in the Minkowski-space effective action can be matched to those in the Euclidean-space effective 
action using a trivial analytic continuation, E -> iS. A thorough discussion is contained in [4]. 
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in the effective action are constrained by the underlying symmetries: gauge invariance, Lorentz 
invariance, "rf, and 3T. The leading operator is the magnetic moment term, N G^ V F^ V N, and 
the corresponding effective Hamiltonian is given by: Hjp = —}X n K-E, where the matrices K are 
generators of boosts in the spin- 1/2 representation of the Lorentz group. Operators at the next 
order contain the s-wave and d- wave couplings of the neutron to two photons, N N F^ V F^ V , and 
i(NYiid v N — d v N y^N)F^ p F p v \ While the latter term would naively be supressed by a power of 
the neutron mass, a field redefinition is required to arrive at a canonically normalized kinetic term, 
and both operators end up being of the same order. A linear combination of these two operators 
gives rise to the second-order effective Hamiltonian: = — ^gceE 2 . 

Both effective interactions, and Hjp, give rise to electric dipole moments (EDMs). The 
magnetic moment interaction generates a motional EDM [5], d^ 1 ' = \l n G x v, where v is a small 
neutron velocity. The interaction energy of the dipole and the electric field is merely — rfW E = 
—Jl • B, where B = v x E is the magnetic field in the neutron's rest frame. The induced EDM, 
(ft 2 ) = — oleE, is proportional to the strength of the applied field with a coefficient that is the 
electric polarizability. 

The electric polarizability lowers the neutron's energy. For a neutron at rest, the energy shift is 
merely AE = — ^cXeE 2 - Less obvious, however, is the effect of the magnetic moment. The leading- 
order contribution clearly vanishes for a neutron at rest. One must treat the magnetic moment 
operator to second order to account for all terms in the neutron energy at 6{E 1 '). For a small 
velocity v, the second order contribution has the form 

oc^l„(vxE)- — j-— -Hn(vxE), (2.1) 
0- ^M n v 2 

and survives the v — > limit. In this way, one avoids the neutron pole. Beyond this schematic 
discussion, one can make the result exact by calculating the neutron propagator to all orders in 

(2) 

and H„ . The result is a neutron propagator with an energy shift of the form 




(2.2) 



As a consequence, studying the electric -field dependence of unpolarized neutron correlation func- 
tions is not enough to determine the electric polarizability without knowledge of the magnetic 
moment. 

In order to determine Cfe from two-point functions, one must perform the analogue of a Born 
subtraction. For the case of a magnetic field, the magnetic moment can be isolated by considering 
spin-projected correlation functions. One way to access magnetic moments in an electric field is, by 
analogy, to use boost-projected correlation functions. We use this terminology to refer to the spinor 
structure; the neutron remains at rest throughout. With projection matrices, £P± = i(l iA^), we 
observe that the boost-projected two-point functions have the form 

Tr[^±G(0] =Z(l±^ n E)e- it{M " +AE \ (2.3) 

for the electric field E = Ez, with the energy shift AE given in Eq. (2.2). Thus a simultaneous mea- 
surement of both boost-projected correlation functions will allow one to determine the magnetic 
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moment and electric polarizability. 2 

While we have focused our discussion on the neutron, a similar analysis is also possible for 
the proton. One might be worried that proton observables would be inaccessible in background 
electric fields because energy is not a good quantum number. This is not a fundamental setback, 
however, as one need not work with states of good energy. We suggested the approach to handle 
charged particles in electric fields using single-particle effective actions [6]. One must sum the 
Born-level couplings to the total charge to arrive at the expected form of the charged particle cor- 
relation function. For the proton in an electric field, there is an additional Born coupling to the 
anomalous magnetic moment, and finally a non-Born term (which is the electric polarizability). 
Nonetheless, the utilization of boost-projected proton correlation functions allows one to access 
the anomalous magnetic moment and electric polarizability. The correlation functions do not have 
simple exponential behavior in time, but can be determined in a model-independent fashion. 

3. Nucleon in Electric Fields on a Lattice 

We now focus on our Euclidean space lattice calculations in background electric fields. To im- 
plement a background electric field, we can choose either a compact or non-compact formulation. 
We utilize the former in order to arrive at uniform fields [7, 8]. Our compact U{\) gauge field is 
implemented with links 3 



Notice the links are unitary as mandated by ^ invariance, i.e. the parameter $ is real valued. With 
qS = 27in/LT, the electric field through each elementary plaquette of the lattice is uniformly S 
in the x$ -direction. Some time ago, we investigated the effects of electric field gradients by using 
non-quantized fields [9]. While energy shifts to neutral particle correlators were found to be on the 
percent level, this is undesirable because the expected energy shifts due to polarizabilities are the 
same size. 

To demonstrate our method of computing nucleon magnetic moments and electric polariz- 
abilities, we used an ensemble of anisotropic clover lattices generated by the Hadron Spectrum 
Collaboration. Table 1 summarizes the lattices used in our computation. Because these lattices 
were designed for the investigation of excited-state hadrons, we utilize two-state fits in order to 
stabilize the extraction of ground-state parameters. We provide a survey of our results [2] focusing 
on three cases: unpolarized neutron, boost-projected neutron, and boost-projected proton. 

For unpolarized neutrons, we can only determine a combination of magnetic moment and elec- 
tric polarizability, Eq. (2.2). To do this, we measure the unpolarized neutron correlation function 
for several values of the external electric field. Typical such measurements are shown in Fig. 1. By 
measuring the neutron energy as a function of the applied field S, we can then extract the coefficient 

2 As the lattice correlation functions are determined in Euclidean space, it is useful to quote the Euclidean version 
of Eq. (2.3): Tr [^±G{t)} = Z( 1 ± \l n S ) exp[-;(M„ + AE)}, where t is now the Euclidean time, &>± are the "boost" 
generators in Euclidean space, and AE = +j ( Ote — j^r ) <=' 2 - 



Currently our simulations are restricted to post-multiplied background fields. As a consequence, physical predic- 
tions can only be made for the isovector nucleon magnetic moment, although the isovector polarizabilities are expected 
to be fairly insensitive to the sea quark charges. 



'n(x) = exp(-iqS'x48^ :3 )exp(iqS , Tx 3 d^ A 8x 4 ,T-i)- 



(3.1) 
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Table 1: Lattices, propagator inversions, and background field strengths used. Further details of the lattice 
action can be found in Refs. [10, 11] with N s — » 20 for this work. The electric field strength is listed in terms 
of the integer n appearing in the quantization condition. 
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Figure 1: Effective mass plots for unpolarized neutron correlators, Tr[G„(f)]. For two values of the electric 
field (corresponding to n = 1 and n = 2), we show the effective mass of the neutron correlator along with 
the effective mass of our two-state fit to the data with uncertainty band. The extracted ground-state energy 
is shown as a (red) flat band. 

of the <f 2 term. This is not simply the polarizability, rather the combination s/ n = | {ole — jf^j- 
Using the unpolarized neutron data, we find si n = 1 .3(9) ( 1 ) ( 1 ) x 10~ 4 fm 3 . 

On the other hand, with the boost-projected correlation functions, we can separate out the neu- 
tron magnetic moment and electric polarizability. This is accomplished by measuring Tr[£P±G n (t)] 
correlators for various values of the applied field (the strength of which is indexed by n), and 
then performing simultaneous fits for each n. Such fits allow one to determine the energy shift, 
AE = s^nS 1 + . . ., and the amplitudes Z(l ± \i n S\ From this information, one can determine the 
coefficient stf n of the £ 2 term in the energy shift. We find stf n = 1.3(7)(2)(1) x 10~ 4 f m 3 , which is 
consistent with our analysis of the unpolarized correlators. Unlike that case, however, we can go 
further and determine the magnetic moment, \i n = —1.63 (10) (4) (5) [/iy], and electric polarizabil- 
ity, a' E = 3. 3(1. 5) (2) (3) x 10~ 4 fm 3 . It should be emphasized that these are "connected" values. 

Finally we can perform a similar analysis for the boost-projected proton correlation functions. 
The time-dependence of the proton correlation functions is considerably different than that of the 
neutron because of the Born couplings to the total charge. While there is no simple exponential 
falloff at long times, there is a model-independent prediction for the behavior of the correlator. Us- 
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Figure 2: Effective mass plots for boost-projected neutron correlators, Tr\3P±G n (t)\. For a single value 
of the electric field (corresponding to n — 1), we show the two effective masses of boost-projected neutron 
correlators. For these two plots, we perform a simultaneous two-state fit in order to extract the neutron 
magnetic moment and electric polarizability. Results of the fit are also shown on the effective mass plots. 
The extracted ground-state energy is also shown as a (red) fiat band. 



P:n=2 



P: n=2 




Figure 3: Effective mass plots for boost-projected proton correlators, Tr[&±G n (t)}. For a single value of 
the electric field (corresponding to n = 2), we plot the two effective masses of the boost-projected correlators. 
The proton correlator is not expected to have a simple exponential behavior at long times. The two effective 
masses of our simultaneous fit to both boost-projected correlators are shown, as well as just the ground-state 
contribution (gray) and ground-state rest energy (red). 

ing the form predicted for a spin- 1/2 particle, we can extract the magnetic moment and electric po- 
larizability by performing simultaneous fits to both boost-projected proton correlators. In Fig. 3, we 
display the behavior of the boost-projected proton correlators calculated in a particular electric field 
strength. From the fits over the various field strengths, we determine }X p = 2.63(13)(1)(4) [jUjv], and 
Off = 2.4(1. 9) (3) (2) x 10~ 4 fm. These are again only "connected" contributions. 

4. Outlook 

We have outlined our progress in the computation of hadronic electromagnetic properties us- 
ing the background field method. In background electric fields, we showed that nucleon correlation 
functions depend on the magnetic moment and electric polarizability. Unpolarized neutron correla- 
tors do not allow one access to the electric polarizability, only the combination AE in Eq. (2.2). By 
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contrast, boost-projected correlation functions allow separate determination of magnetic moments 
and electric polarizabilities. Using an ensemble of anisotropic gauge configurations (courtesy of the 
Hadron Spectrum Collaboration), we demonstrated our method successfully for both the neutron 
and proton. The proton had yet to be treated in background electric fields. 

There are many refinement possible to our computation. The usual caveats must be issued: 
we have computed quantities at only one value of the pion mass, and one lattice spacing. Due to 
computational restrictions, the background field was not included in the gauge field generation. 
Additionally there are non-standard finite volume effects that one must worry about [12, 13, 14]. 
For example, the uniform electric field does not, strictly speaking, lead to a periodic lattice action. 
The action is only periodic up to a gauge transformation. Charged particle correlation functions 
experience a gauge defect at the edge of the lattice. This, however, is a calculable effect. One 
can determine the effect of the defect on the two-point function by using a compact single-particle 
effective action. Another finite volume oddity concerns the non-trivial holonomy of the background 
field. New interactions are generated from virtual particles wrapping around the lattice. These 
should be exponentially small, moreover their effects can be addressed in an effective field theory 
framework. Work is currently underway to address these issues. 
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